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1. Introduction 

Although lattice calculations have given us much insight into the non-perturbative regime of QCD, 
a thorough understanding of infrared phenomena like confinement and chiral symmetry breaking 
has not been achieved yet by lattice calculations as even the largest lattice calculations presently 
available do not yet probe sufficiently deep into the infrared. Hence, in recent years there have 
been many activities devoted to the non-perturbative studies of the infrared sector of continuum 
QCD. Among these are functional renormalization group (FRG) equations [[I], and the variational 
approach to the Hamiltonian formulation of QCD [||, Q|. Each approach has its own advantages 
and drawbacks and by combining these approaches one can expect to gain new insights into the 
theory, in particular into the non-perturbative regime. 

My talk is devoted to the application of FRG flows to the Hamiltonian formulation of Yang-Mills 
theory as put forward in First, to motivate the FRG studies, I will briefly summarize some 
essential results of the variational Hamiltonian approach. Then I will give a short general introduc- 
tion to FRG equations and present the new flow equations for Hamiltonian Yang-Mills theory in 
Coulomb gauge together with their solutions for the gluon and ghost propagators [[Tj]. 



2. Summary of the variational approach to Hamiltonian Yang-Mills theory 

In the variational approach to Yang-Mills theory in Coulomb gauge developed in our group one 
makes the following ansatz for the vacuum wave functional [Q] : 

Y = , 1 e*p \-\ [ AcoA] , (2.1) 

^Det(-Dd) V 2J J 



^Det(-Dd) \ 2J J 

where Det(— Dd) is the Faddeev-Popov determinant and ft) is a variational kernel, which is related 
to the static gluon propagator by 

(AA) = (2(a)- 1 (2.2) 

and therefore has the meaning of the gluon energy. Minimizing the vacuum energy density with 
respect to co(p) [Q] one finds the gluon propagator shown in Fig. |] together with lattice data 
[Eh. Both, the lattice data and the continuum solution can be fitted with Gribov's formula 



M 4 

with M ~ 860 MeV . At large momenta co(p) ~ p, in accordance with asymptotic freedom, whereas 
at small momenta co(p) ~ l/p, and this divergence shows the absence of gluons from the physical 
spectrum, also called gluon confinement. Fig. ^| shows the ghost form factor d(p), defined by 

((-Dd)- l ) = 4£ 1 (2.4) 
P 

obtained in the variational approach and on the lattice. In the variational approach one can show 
that the infrared exponents of the ghost and gluon propagators, 

G>(p^0)~l/ p a , d(p^0)~l/pP, (2.5) 
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Figure 2: Ghost form factor from the variational 
Figure 1: Gluon propagator from the variational approach compared t0 lattice data 

approach compared to lattice data 

are related by a sum rule under the assumption of a trivial scaling of the ghost-gluon vertex pfl, 

a = 2j3 - 1 . (2.6) 

Furthermore, assuming an infrared scaling behaviour of the ghost form factor, i.e. j8 > 0, one finds 
two "critical" solutions 

i) a = 0.6 , jS = 0.8 it) a = 1 , /3 = 1 . (2.7) 

The latter results in a linearly rising quark potential shown in Fig. ||. However, sub-critical solu- 
tions d~ l (0) ^ do also exist [H]. Previous lattice calculations in Coulomb gauge seemed to yield 
j8 = \. However, these calculations were carried out on rather small lattices. As one increases the 
lattice size, j8 seems to increase as well and recent lattice calculations |Q] are at least compatible 
with j8 = 1 . In Coulomb gauge the inverse ghost form factor has been shown to represent the di- 
electric function of the Yang-Mills vacuum e(p) = d~ l (p), and the so-called horizon condition 
d~ l (0) = implies then that the Yang-Mills vacuum is a perfect dual color superconductor. 
The advantage of the Hamiltonian formulation is its close connection to physics. In the variational 
approach one makes an ansatz for the unknown vacuum wave functional which encodes all the 
physics. This ansatz can be systematically improved towards the full theory. The price to pay is 
the apparent loss of renormalization group invariance which already complicates the discussion of 
perturbative renormalization, but also the discussion of scaling laws in the infrared. 

3. FRG flow equation 

Renormalization group invariance is naturally built-in in the functional renormalization group ap- 
proach to the Hamilton formulation of Yang-Mills theory which hence has the advantage of com- 
bining renormalization group invariance with the physical Hamiltonian picture. 
The starting point of the FRG flow equation approach [Q] is the (renormalized) generating func- 
tional of Green's functions 

Z[j] = J 9(pe- s ^ + ^ , (3.1) 
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where (p and j denote collectively all fields and sources. In the FRG approach Z[j] is infrared 
regulated by adding the regulator term 



r 1 1 1 f 

A $k[(p] = -(p-R k -<p = -J ^ Rk< ^ 



,-S[q>]-AS k [q>}+j-(p = Mj] 



to the classical action, 



The regulator function R k {p) is an effective momentum dependent mass with the properties 

Urn R k {p) > , lim R k (p) = , 

p/k->0 k/p—>0 



(3.2) 



(3.3) 



(3.4) 



which ensures that R k (p) suppresses propagation of modes with p < k while those with p > k are 
unaffected and the full theory at hand is recovered as the cut-off scale k is pushed to zero. 



By taking the derivative of Eq. (3.3) with re 



spect to the (dimensionless) momentum cut-off 
t = In one obtains the flow equation for 
the generating functional of connected Green's 
functions 

3BrM \8W k .8W k l T .5 2 Wt 

(3.5) 

(where the dot denotes d t ). In practice, it is 
usually more convenient to perform first the 
Legendre transform from the source j to the 
classical field 




Figure 3: Static quark potential 







8W k [j] 
Sj 



(3.6) 



resulting in the effective action 

r k [(j>) 



-w k [j]+j-(!>) j=m --<!>-R k -<l> , 



(3.7) 



where j k [(j)] is given by solving Eq. (3^) for j. By taking the derivative of Eq. ( p?7| ) w.r.t. t = \nk/ko 



and using Eq. (3.5) one derives Wetterich's flow equation for 1^, 



where 



d t r k [(j)) 



i 



-Tr- 



1 



5 n r k [$) 



Rk, 



(3.8) 



(3.9) 



501 .. . 8<$> n 

are the one-particle irreducible «-point functions (proper vertices). The generic structure of the 



flow equation (3.8) is independent of the details of the underlying theory, i.e., of the explicit form 



of the action S[0], but is a mere consequence of the form of the regulator term (3.2), i.e., that it is 
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quadratic in the field. By taking functional derivatives of Eq. (3.8) one obtains the flow equations 

ads 

(3.10) 



for the (inverse) propagators and proper vertices. For the two-point function this equation reads 

1 



a r (2) _ I Tr b ; 1 J _ r (4) , 

2 n Z) +R k 



rg ' i£>-hk-2, 



where all cyclic indices (summed over in the trace) have been suppressed. 

4. Hamiltonian flow 

In the Hamiltonian approach the generating functional of static correlation functions reads 

Z[j] = • (p)\y) = J @(p\y[(p}\ 2 exp(j ■ <p) , (4.1) 

where (<p|yf) = y/[(p] is the vacuum wave functional. With the identification 

|i/%]| 2 = exp(-%]), (4.2) 

defining the classical action S[tp] in terms of the vacuum wave functional y/"[<p], Eq. ( p~H ) has 
precisely the standard form of a generating functional except that the functional integral extends 
over time independent fields. Here we are interested in the Hamiltonian flow of Yang-Mills theory 
in Coulomb gauge whose generating functional reads 

Z[J] = J ^ADet(-D5)|v/[A]| 2 exp(7-A) , (4.3) 

where the integration is over transversal gauge fields A and the Coulomb gauge condition has been 
implemented by the usual Faddeev-Popov method. Representing the Faddev-Popov determinant in 
the standard fashion by ghost fields, c,c, 

Det(-Dd) = J »c^ /f( - M)f (4.4) 

the underlying action, cf. Eq. (fO|), reads 

S[(p] = -ln\y[A]\+ Jc(-Dd)c . (4.5) 

The general flow equation ( |3.10| ) still holds provided that (f> is interpreted as the superfield (f> = 
(A,c,c). The FRG flow equations for the gluon and ghost propagators are diagrammatically given 
in Fig. | 

5. Approximation schemes and numerical solution 

The FRG flow equations embody an infinite tower of coupled equations for the flow of the prop- 
agators and the proper vertices. These equations have to be truncated to get a closed system. We 
shall use the following truncation: we only keep the gluon and ghost propagators, to wit 

„2 



(3) 

In addition, we keep the ghost-gluon vertex i k i dc , which we assume to be bare, i.e., we do not solve 
its FRG flow equation. The above Uuncation removes the tadpole diagrams from Fig. |[ Moreover, 
we shall assume infrared ghost dominance and discard gluon loops. Then the flow equations of the 
ghost and gluon propagator reduce to the ones shown in Fig. |5[ 
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Figure 4: Flow equations of the propagators. The spiral and dotted lines with black circles denote the 
regularized gluon and ghost propagators at cutoff momentum k, respectively. White circles stand for proper 
vertices at cutoff k, a regulator insertion R k is represented by a square with a cross. 
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Figure 5: Truncated propagator flow equations. The bare vertices at A: = A are symbolized by small dots. 



The flow equations in Fig. |5] are solved numer- 
ically using the regulators 



R A ,k(p) 

nip) = 



2prk(p) , R c ,k(p) = p 2 n(p) 



exp 



n 2 



El 

k 2 



(5.2) 



1e-02~- _ 

1 e+01 „ 
p 1e+04^ 



1e+04 1 e+02 



1e+00 1e-02 1e-04 



and the perturbative initial conditions at the 
large momentum scale k = A, 

K 

d A (p) = d A = const . , (0 A (p)= P + a . (5.3) 

Figure 6: Flow of the ghost form factor. 

With these initial conditions, the flow equations 

for the ghost and gluon propagators are solved under the constraint of infrared scaling for the ghost 
form factor. The resulting full flow of the ghost dressing function is shown in Fig. |6[ As the IR 
cut-off momentum k is decreased, the ghost form factor dk(p) (constant at k = A) builds up infrared 
strength and the final solution at k = k m i„ is shown in Fig. [7] together with the one for the gluon 
energy (Ou min {p)- ^ i s seen tnat th e IR exponents, i.e., the slopes of the curves dk mjn (p),(Oic mjn (p) 
do not change as the minimal cut-off & m! „ is lowered. Let us stress that we have assumed infrared 
scaling of the ghost form factor but not the horizon condition. The latter was obtained from the 
integration of the flow equation but not put in by hand (the same is also true for the infrared 
analysis of the Dyson-Schwinger equations following from the variational Hamiltonian approach, 
i.e., assuming scaling the DSEs yield the horizon condition). The infrared exponents extracted 
from the numerical solutions of the flow equations are 



a = 0.28, [5=0.64 



(5.4) 



They satisfy the sum rule <\2.6\) resulting from the DSE obtained in the variational approach but are 



smaller than the ones of the DSE, see Eq. < \2.T\ ). Moreover, the present approach allows to prove 
the uniqueness of the sum rule 



3) |Jl|l, analogously to the proof in Landau gauge [jlOJ. 
Replacing the propagators with running cut-off momentum scale k under the loop integrals of the 
flow equation by the propagators of the full theory, 



d k (p) ->■ d k=Q (p) , co k (p) ->■ (D k = (p) , 



(5.5) 



6 



Hamiltonian Flow of Yang-Mills Theory in Coulomb Gauge 



Hugo Reinhardt 



10 




-7 -5 -3-11 3 1U -7 -5 -3 -1 1 3 

10 10 10 10 10 10 10 10 10 10 10 10" 

p [GeV] p [GeV] 



Figure 7: Inverse gluon propagator (O and ghost form factor d at three minimal cutoff values £,„,„. 

amounts to taking into account the tadpole diagrams [jl|]. Then the flow equations can be ana- 
lytically integrated and turn precisely into the DSEs obtained in the variational approach to the 
Hamiltonian formulation of Yang-Mills theory [Q], with explicit UV regularization by subtraction. 
This establishes the connection between these two approaches and highlights the inclusion of a 
consistent UV renormalization procedure in the present approach. 

The above results encourage further studies, which includes the flow of the potential between static 
color sources as well as dynamic quarks. 
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